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Abstract. For a static, perfect fluid sphere with a general equation of state, 
we obtain the relativistic equation of hydrostatic equilibrium, namely the 
Tolman-Oppenheimer-Volkov equation, as the thermodynamical equilibrium in 
the microcanonical, as well as the canonical, ensemble. We find that the 
stability condition determined by the second variation of entropy coincides with 
the dynamical stability condition derived by variations to first order in the 
dynamical Einstein's equations. Thus, we show the equivalence of microcanonical 
thermodynamical stability with linear dynamical stability for a static, spherically 
symmetric field in General Relativity. We calculate the Newtonian limit and find 
the interesting property, that the microcanonical ensemble in General Relativity 
transforms to the canonical ensemble for non-relativistic dust particles. Finally, 
for specific kinds of systems, we study the effect of the cosmological constant to 
the microcanonical thermodynamical stability of fluid spheres. 
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1. Introduction 

The deep connection between gravity and thermodynamics is being extensively studied 
for several decades now, starting from the four laws of black hole mechanics [TJ [2] and 
Hawking radiation |3] . These laws were derived by Einstein's equations. More recently, 
by reversing the logic around, many attempts have been made to derive Einstein's 
equations from thermodynamics [U [5] [6] and thus realize gravity as an emergent 
phenomenon, rooted in a deeper, underlying and yet unknown, microscopic structure. 
All these attempts were based on thermodynamics on an horizon or on some other 
boundary, making always use of some form of holographic principle, pointing in this 
way to a quantum gravity theory. However, we think it is crucial, before introducing 
some sort of exotic or quantum principle, to understand at purely classical level [7], as 
deeply as possible, the thermodynamics of self-gravitating gas. In this way, it will be 
possible to figure out which peculiar properties of gravitational thermodynamics are 
of quantum nature, rooted in a deeper theory, and which are just built-in features of 
classical General Relativity. In this spirit, we will study the classical thermodynamics 
of self-gravitating gas in General Relativity with no reference to the thermodynamics 
on an horizon or a screen. 

In fact, the attempt to derive Einstein's equations from thermodynamics of 
self-gravitating gas dates back to 1965 in a work of Cocke [8 1. He proposed a 
maximum entropy principle for a static, spherically symmetric, perfect fluid to derive 
the relativistic equation of hydrostatic equilibrium, that is widely known as Tolman- 
Oppenheimer-Volkof (TOV) equation. In his formulation the fluid was in adiabatic 
motion so that the total entropy was constant and the initial data constraint equation 
together with the radial-radial component of Einstein's equations were used. Later, 
in 1981, Sorkin, Wald and Zhang (SWZ) [9] developed a different entropy principle 
for radiation, which did not need the adiabatic condition and used only the constraint 
equation. Very recently, Gao [10] generalized SWZ's formulation for an arbitrary 
perfect fluid. However, in Gao's 10] formulation is not clear in which thermodynamic 
ensemble the calculation is performed. It is crucial to specify the ensemble, because in 
Gravity the thermodynamic ensembles are not equivalent [11] 112] regarding stability 
properties. We will sharpen his argument, by deriving the TOV equation from 
thermodynamics, working clearly both in the microcanonical as well as the canonical 
ensemble and in a hopefully, more straightforward and simply way. 

Most importantly, we show that gravitation is equivalent to thermodynamics in 
the microcanonical ensemble, not only regarding the equilibria (i.e. TOV equation), 
but even at the level of stability. Indeed, we find that linear dynamical stability, defined 
by Einstein's equations, coincides with microcanonical, thermodynamical stability, 
defined by the second variation of entropy, for a static, perfect fluid sphere with a 
general equation of state. In addition, we show that the microcanonical ensemble in 
General Relativity, becomes the canonical ensemble in the limit of non-relativistic, 
dust matter. 

Finally, we study the effect of the cosmological constant on the stability of fluid 
spheres. This study can be considered as the relativistic generalization of previous 
studies of self-gravitating gas in the Newtonian limit [HI [X5l US]. We perform the 
analysis for matter and radiation with a linear equation of state with the assumptions 
of constant entropy per baryon and chemical composition throughout the whole 
fluid. A general result is that an increase in the cosmological constant tends to 
stabilize the system. In addition, the critical minimum radius, down to which there 
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exist equilibrium configurations, 'hits' the cosmological horizon at some value of the 
cosmological constant. 

Our work is organized as follows. In section [5] we derive the condition for linear 
dynamical stability from the Einstein's equations. In section [3] we prove that TOV 
equation can be derived from the microcanonical ensemble, i.e. by the first variation of 
entropy with fixed mass and number of particles, and we discuss the canonical ensemble 
as well. In section |4] we calculate the second variation of entropy and prove that the 
condition for microcanonical, thermodynamical stability is equivalent to the condition 
for linear dynamical stability calculated in section [5] In section [5] we calculate the 
Newtonian limit of the condition for microcanonical thermodynamical stability, found 
in section |4] Finally, in section [6] we study the effect of the cosmological constant to 
the stability of fluid spheres with specific equations of state. 



2. Linear dynamical stability 

In case of a perfect fluid that is static and spherically symmetric, the Einstein equations 
reduce to the equation known as Tolman-Oppenheimer-Volkoff (TOV) [T7J [TH] (see 
[T5] for a study on TOV with a cosmological constant). TOV is the equation that 
describes hydrostatic equilibria. In |Appcndix A| we review how this equation is derived 
from Einstein's equations. In this section we want to find the condition for linear 
dynamical stability, i.e. for stability under radial variations, to first order, about the 
equilibrium. This problem has been studied by Yabushita [20) . In |Appcndix A| we 
review analytically the calculations of Yabushita, including a cosmological constant, 
as well. Our final purpose is to compare dynamical stability with thermodynamical 
stability, which we will finally do in section [4] For this we need the expression for 
linear dynamical stability of a general equation of state which we calculate in the 
followings. 

The metric for a spherically symmetric system can be written in the form [17] 

ds 2 = e v c 2 dt 2 - e x dr 2 ~ r 2 (d0 2 + sin 2 9d 2 <p) (1) 
with v = f(r,t) and A = X(r,t). For a perfect fluid, the energy momentum tensor is 

TZ = {p + P c 2 )g av ^^~p5t (2) 
ds as 

which at the equilibrium becomes just 

T u = diag(/? e C 2 , -p e , -p e , -p e ) (3) 

where the suffix e denotes quantities at the equilibrium. At the equilibrium the metric 
is static and the Einstein's equations reduce to just two equations (see | Appendix A| : 

, ,p e . ( Gmir) „ » e A 2 \ / 2Grn(r) A 2 \ _1 fA . 

m! = iirr 2 p e (5) 

The first one is the TOV equation that defines hydrostatic equilibria, while the second 
implies that m(r) is the total mass-energy inside radius r [21] . The equilibria can be 
determined by these two equations, provided an equation of state p = p(p) and initial 
conditions p(0), m(0) are given. 

Now, let us apply small perturbations about the equilibria (see |Appendix A[ ) as 
in equation (|A.14[) . Let 

f(r,t) = 6m(r,t) (6) 
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with 

Then, assuming a perturbation 

6m ~ e at (8) 
with Sm(r = 0) = one obtains from Einstein's equations, to first order, that: 

¥ + ^-(5 P + s P c 2 ) + S -^( Pe + Pe c 2 ) = e -^r° 2 f (9) 

This is equation (2.13) in Yabushita [20] and is proved in |Appendix A with a 
cosmological constant included in the calculations. It is clear from equation (j8]) that if 
a 2 < the system will pulsate about the equilibrium and therefore it is stable. Given 
an equation of state, all variations in equation ((9|) can be expressed with respect to 
/, so that (0) becomes a differential equation with respect to / and a differential 
operator may be defined on the left hand side as we will see. The stability is therefore 
determined by the sign of the eigenvalues of this operator. If it has any positive 
eigenvalue the system is unstable. 

To first order in perturbations (|A. 14|) . equation (|A.2I) becomes 



Substituting 



c 4 ^ \ rc 2 3c 2 ^ A / 1 r 



2Gm 8ttG „\ 1 2G 



rc 2 3c 2 

solving with respect to Sv' and substituting it together with z/, from equation (|A.7|) . 
in equation ([9]) we get: 

x' P e ' tx , x 2n , , , 2 A (, 2Gm 8nG 2 X 

dp ; j \°P + °P C ) + CPe + PeC ) 1 5 VirP Kr 

Pe + PeC V rc ^ c 

This is an important equation for our purposes, since it depends only on 5p and 
/ = 8m. This enables us to compare it with the corresponding expression deduced 
from thermodynamical stability and in this way determine the necessary condition for 
the equation of state, so that the two kinds of stability to be equivalent. This will be 
done in section |U 

Let us assume a general equation of state of the type 

p = p(p) 



Let 



so that 



and 



, , 1 dp , . 

^ = ~ 2 ip M 

, gc 2 df 

p' = g(r)c 2 p and dp = g(r)c 2 Sp = - — -— (12) 

47rr z dr 

r / c 2 / d 2 f , / , 2g\ df\ 
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Substituting equations (fT2|) . (flU)) and (JT} into (|TU)) we finally get: 



x 1 



2G 



x 1 



re" 
2Gm 



8ttG 
Ite 2 

SttG 



d 2 / 
dr 2 

PAT 2 



2 (g + l)cV 

J" Pe + PeC 2 



4ttG 



rr 



3c 2 



(p e + PeC 2 ) 

2gc 2 p e ' 



'1 

dr 



AmG . 2 , 
H 5~ r ( Pe + ^ eC ) 



(14) 



p e + p e c z 

Thus, linear dynamical stability is determined by the sign of the eigenvalues of the 
operator 



4ttG 



dr 2 



(ff + l)cW , 4ttG 
5 1 —r 

p e + p e C 2 



2Gm 8ttG 



3c 2 



PAT 



-1 



H 2~(Pe + /°e c ) r 1 



2Gm 8ttG 



3c 2 



{Pe + p e C 2 ) 

2gc 2 p e ' 



d 
dr 



Pe + p e C 2 



(15) 



with g = dp/dpc 2 . If there exists any positive eigenvalue the equilibrium is unstable. 
In a series of equilibria, a turning point of stability happens at the equilibrium for 
which: 



Lf = 



(16) 



3. General relativistic equation of hydrostatic equilibrium from 
thermodynamics 

As we have stated in the introduction, Gao [10] generalized SWZ's [9] formulation 
and managed to derive TOV equation by extremizing the entropy, for an arbitrary 
perfect fluid, but with no reference to any specific thermodynamic ensemble. We 
know that specifying the ensemble is very important in gravity, because they are not 
equivalent. In this section we will derive the TOV equation from thermodynamics, 
working clearly in the microcanonical ensemble and performing the calculation in 
a different, straightforward way. We will include in the calculation a cosmological 
constant, as well. Next, we will discuss the canonical ensemble, too. 

Following Gao [TU] and SWZ [S], let us state the main assumptions: 

• We assume a perfect fluid for which the first law of thermodynamics holds: 

c 2 

ds = —dp - —dn (17) 

where s is the entropy per unit volume, T is the temperature, p is the mass 
density, related with the energy density by e = pc 2 , fi is the chemical potential 
and n is the baryon number density. This equation implies that p and n can be 
considered as two independent variables and s — s(p,n). 

• The Gibbs-Duhem relation, which is a general expression for a perfect fluid, holds: 

Ts = pc 2 + p — pn (18) 
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• The system is spherically symmetric and the entropy maxima necessarily 
correspond to static configurations!]]. Therefore, the initial value constraint 
equation [3] for time symmetric data holds: 

Wr= 1 -^( p + Pa ) (19) 

where ^ R is the intrinsic curvature, p is the matter density and /?a is the density 
corresponding to a (positive or negative) cosmological constant A 

Ac 2 , s 

' A = 8^G (20) 
In terms of the radial-radial component of the metric g rr we have 

m *=^i{^-9rr 1 )} (2D 



Using equation (|2T|) to integrate ([T9]) . we get: 

2Gm(r) 8ttG , 



9rr = [I - - ^PA^j (22) 

with 

m(r) = j p(r)4xr 2 dr (23) 
Jo 

These are our assumptions. Let R denote the radius at which the pressure vanishes 
[18], so that it defines the edge of the fluid sphere. Using spherical symmetry and 
equation (|22|) in the proper volume element, the entropy can be written as: 

f R , ( 2Gm(r) 8nG ,\~* . 2, 
S = J s(r) ( 1 - ~^ PAT ) 4nr dr (24) 

Similarly, the baryon number is given by 

f R 2Gm(r) 8ttG 2 \~' 2 , 

A^ - / n(r) 1 ^ - -^PAr 2 4irr 2 dr (25) 

J \ rc Ac J 

In contrast, to obtain the total mass-energy M, including that of the gravitational 
field, one should not integrate on the proper volume but on the normal volume (see 
p. 302 in Weinberg [21]) as in equation ([23]) 

nR 

M = m{R) = / p{r)ATir 2 dr (26) 
Jo 

We want to find the extrema of entropy under the constraints of fixed energy and 
baryon number, so that we work in the microcanonical ensemble. We will use the 
method of Lagrange multipliers. Let /3o, a be two, yet undetermined, Lagrange 
multipliers. The requirement of entropy extrema under the constraints of constant 
energy and baryon number is given by 

SS - Pqc 2 6M + aSN = <=> 
5s(g rr ) 2 r dr H — j sSm(g rr ) 2 rdr — (3qc / Spr dr 



c 2 



G r R 



aSn(g rr ) 2 r 2 dr H — - / an5m,(g rr ) 2 rdr — (27) 
I Note that this is proved by SWZ [9] for the asymptotically flat case. 
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Ss=pp-^8n (28) 

in the first term and Gibbs-Duhem relation (jl8p to substitute s in the second term, 
equation (f2~T|) becomes 



R 1 w ^ 2J G f R pc 2 + Pl ^ 



c / 7^ S P(9rr) 2 r dr + — j — 5m(g rr ) 2 rdr - /3 c 2 \ 5pr 2 dr 

,-R 



f G f 3 

+ J (a-^)6n(g rr )K 2 dr+^j (a - n5m(g rr )§ rdr = (29) 

Since p and n are independent variables and Sm ~ <5p, the fourth term of the last 
equation ([29]) implies 

a = £ (30) 

This is an expected result, since the temperature and chemical potential should suffer 
the same red-shift. Using this value for a and equation (|23l) we get 



^;Sp{g rr ) 2 r 2 dr + ^ f ^-^ P ' 



c / 7^$p{9rr) 2 r dr + — j — - — {g rr ) 2 r[l S p(r)4nr df ) dr 

-R 



-/3 c 2 / 5pr z dr = 
Jo 

: <MM(cV) • + ^ jf* ^+^) ( grr (r))l r -df - A, | = (31) 

where, going from the first line to the second, we have interchanged in the double 
integral the integration variables, changing the limits of integration accordingly. In 
order for this equation to hold for any variation 5p, it should be 



1 a 4ttG f R p(?)c 2 +p{r) 



3 



^ = ^{9rrY + —J T (f) \9rr{f)Yfdf (32) 

which for To = l//?o an d r — R becomes 

To = T(R)(g rr (R))-i or T = T(R) (l - ^ - ^W? 2 ) ' (33) 

which is the well known Tolman relation |17) . Therefore the Lagrange multiplier 
To = l//?o is the surface temperature at the edge r — R, measured by an observer at 
infinity. 

Let us repeat the proof of a relation given in Gao's Addendum [TU], that we will 
use here, as well. By differentiating the Gibbs-Duhem relation (ITS! and using the first 
law ([!?]) , we get 

dp = sdT + ndp p' = sT' + np' (34) 
Equation ([3171) gives 

p! = aT' (35) 
Combining equations (j34)) . (|35|) and (|T8|) we have 

r = (36) 
1 p + pc z 
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This is the desired equation. Now, differentiating equation we get 
T ( 2Gm(r) 8nG 2 \~? 1 ( 2Gm(r) SttG 



T 2 \ rc 2 3c 2 ^ A ) T \ rc 2 3c 2 ^ A 

I f . „ Gm 8nG 
x — I -AnGpr + — — p\r 

~^0L r + teC( ,> (l - ^2 - - (37) 

which after replacing equation (|36p gives finally the TOV equation 
, ,p , j Gmjr) i A „p 8irG W 2Gm(r) 8^G 2 \ 1 ^ 

Let us discuss the case of the canonical ensemble. Imagine the fluid sphere to be 
bounded by reflecting and insulating walls and let on the outside be a heat reservoir 
with constant temperature To at the outer surface of the walls, as measured by an 
observer at infinity. Since the exchange of heat energy between the system and the 
reservoir happens only at the edge, the temperature that will enter the expression of 
the Helmholtz free energy F is To: 

F = Mc 2 - T S (39) 

The minima of F will give the stable equilibria. Sometimes is useful (and it is 
equivalent) to use the Massieu function J = —/3qF, i.e. 

J = S- /3 Mc 2 (40) 

Now, the maxima of J will give the stable equilibria. Hereafter, we will refer to the 
expression (|40[) as the free energy. Since (3 is constant and is the inverse temperature 
at the edge of the sphere as measured by an observer at infinity, it is trivial to see that 
the extremum, i.e. the point of elimination of the first variation, of the free energy 
with constant baryon number is equivalent to the extremum of entropy with constant 
energy and baryon number 

5 J + aSN = <3> SS - (3 c 2 5M + aSN = (41) 

Thus, the two ensembles, canonical and microcanonical, give the same equilibria, i.e. 
TOV equation, as they should do. But, the second variation of J and S are different. 
So that, the stability properties of the two ensembles are different, as we will see in 
the next section. 



4. Thermodynamical stability 

It is well known that in gravity the thermodynamic ensembles are not equivalent 
[TT1 [T^] . This means that although they describe the same equilibria, the stability 
properties arc different. In general, self- gravitating systems are more unstable in 
the canonical ensemble rather than the microcanonical. This means that in the 
microcanonical ensemble there are more equilibrium configurations rather than in the 
canonical. On the other hand there is dynamical stability, which at least conceptually 
is different. Since thermodynamical stability in the two ensembles is different, the 
dynamical stability has to be equivalent either with only one of them or with no one. In 
the case of Newtonian gravity, for a spherically symmetric system, it has been proved 
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by Chavanis |13j that dynamical stability is equivalent with thermodynamical stability 
in the canonical ensemble. However, in General Relativity, SWZ [5] have proven for 
the specific case of radiation that dynamical stability coincides with thermodynamical 
stability in the microcanonical ensemble. Consequently, it is evident that at least 
in gravity it is not trivial to conclude on the relation between thermodynamical and 
dynamical stability. The question raised and addressed in this section is if dynamical 
stability coincides with microcanonical thermodynamical stability for any equation of 
state of a perfect fluid in General Relativity. 

In the microcanonical ensemble, the stability of a equilibrium is determined by the 
second variation of entropy. If it is negative, it corresponds to a local entropy maxima 
and the equilibrium is stable, or better metastable in our case since the maxima will 
be only local entropy maxima. Of course, the global entropy maxima is given by the 
Bekenstein-Hawking bound at semi-classical level, while at classical level there is no 
global entropy maxima. Let us calculate the second variation of entropy. 

Recall that (equations @, ©, ([2"2"])1 

1 df ( 2Gm{r) 8nG 2 



f = 6m > 6p= 4^dr> 9rr= V-^^-l^^ I (42) 



and 5m(0) — 8m{R) = 0. Also, equation (|3"o]l implies 

ST _ 8p 
T p + pc 2 

so that 

1 \ 1 5p 

T J T p + pc 



(43) 



*U)=-^b (44) 



Using equations (l4"2"j) , (|4"4")l . the second variation of entropy given in (|2"4"|) is: 

fiS= -J Q R dr^^f(g rr )i (T) 

+ I Q R dr^lff r ( 9rr )i (I 2 ) (45) 

+ f R dr& (p + p C i)^P(g rr )l (J 3 ) 

where all quantities (except from variations) correspond to an equilibrium, so that the 
suffix e used in section [5] is suppressed. Integrating by parts integral (I2) we get 

rR 1 - G , ^ f P' 1 o 1 ( a r, Gm SttG \ \ ^ 

^(drr) 2 <—^ + -+3- i-^Gpr+— — p A r )" 1 ' 

lr y p + pc r c \ r z 6 / 

Therefore 

h + h= / dr—{g rr y — ^ + - + 3- 4^G^r + — —p A r g rr } f 2 

Jq lr [p + pc r cr \ c 

which, by use of TOV equation (|A.12I) becomes: 

h + h = ! dr^{g rr y (--^ + -)f ( i<i ! 

J Tr \ p + pc 2 r 1 

Integrating by parts the integral (Ii) we get 
h= I dr^ - \ 2 (grr)* \sp -5p 2p ^ P 2 \f 







t p + pc 2 y p + pc 



R c 2 1 . .3 1 / Gm 8nG , 

dr 7f;— 2 Wrr ) 2 op— -AnGpr + — p\r f 

T p + pc 1 cr \ V s - ' 
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which, by use of TOV equation (|A.12I) in the second line, becomes: 



h = 



dr- 



1 



T p + pc 2 



(g rr ) 2 \ 5p - Sp 



p' + p'c 2 
p + pc 2 



. 47TG. 2 , 
0P^-(P + PC )9r 



(47) 



Substituting equations (j46j) . (j47|) into (|45| we get 



x < <5p' — (5p 



6 2 S 



p'c 2 



dr- 



1 



T p + pc- 



p + pc 2 



(p + pc 2 )g r 



dp 



AttG G 



2p> 



p + pc 2 



f 



(48) 



c rc* 

Comparing this with equation (|10p we see that for the quantity in brackets in the 
above equation to be equal to the left-hand side of (ITUl) . and thus the microcanonical 
thermodynamical stability to be equivalent to linear dynamical stability it must hold: 

p'Sp = 5pp' (49) 

For a general equation of state p = p{p) the above equation holds. Setting g = dp/dpc 2 
it is straightforward to check that equation (|4"5]l becomes: 



dr 2 



6 2 S = 

g r 
AttG 



dr 



1 



TAirr 2 p 
(.9 + l)cV 



:(9rr) 2 f 



p + pc* 

2\ 



pc* 
4nG 

-^rg rr {p- 



d 
dr 



(p + pc )rg. 



2gcV 
p + pc 2 



(50) 



The sign of S 2 S is therefore determined by the sign of the eigenvalues of the operator 



L 



dr 2 



+ - 



g r 
4ttG 



{g + l)c 2 p' 4ttG 



p + pc* 
(p + pc 2 )rg rr 



rg rr (p + pc 



d_ 

dr 



2gc 2 p' 



(51) 



c* \ p + pc* 

This operator is exactly equal to the one in equation (|15p. The equilibrium is stable 
if L has only negative eigenvalues since then it corresponds to an entropy maxima. 
The same was true for L in the dynamical case in section [3J Thus, the condition 
for microcanonical thermodynamical stability is equivalent to the condition for linear 
dynamical stability. 

Let us investigate the canonical ensemble as well. For p = p(p), the second 
variation of the free energy (|40f is: 



S 2 J 



G 



■ df, 



dr—(g + 2)ff(g rr y* 
lr dr 



pc 



(52) 



Calculating the first integral by integrating by parts and substituting it in the above 
equation we get 



6 2 J = 



G 

dr — (q r 
Tr 



2r 



2 3g (Gm 8ttG 

- A grr 1 



2c 2 



Gi 



H — ~g rr 47rG^-r H — 



8ttG 



-PIsJ 



3 

AnGpr 



-PhJ 



(53) 
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Since g rr > (so that no black hole is formed) we see that for a linear equation of 
state 

p = qpc 2 , q = const. 

and for A < it is S 2 J > for any perturbation /. Therefore, in this case, the 
canonical ensemble is completely unstable. 



5. Newtonian limit 



The Newtonian limit of TOV equation is Emden equation and it is taken assuming 
dust particles [301 HI- I n this case, the equation of state is linear with 
kT 



V - 



where m p is the mass of one particle. Defining 



q = 



kT 
m p c 2 



the equation of state can be written as 

p = qpc 2 

and the non-relativistic dust particles limit corresponds to q — >• 0. 

Let us calculate the operator L (f5"Tj) in the Newtonian limit. Usin 



(54) 



(55) 



(56) 



g 9r 



1. 



equations (I55[) . (|56[) in order to have only q and not c in L, and then taking the limit 
q — > 0, we have: 

2 
r 



kT 

m„ 



dr 2 



P_ 

P 



AirGm p 

i — n^ r P(i 



kT 



1) 



d_ 

dr 



■47rG(g + l)rp 



1 



(q+l)p r 



Airpr T 



f k d 


( 1 




G \ 


\ m p dr 


\4irpr 2 


dr J 


+ Tr i] 



(57) 



Therefore, in the Newtonian limit, stability is determined by the eigenvalues of the 
operator 



f k d 
Ln — — 



1 



G 
Tr 2 



(58) 



m p dr \4irpr 2 dr / 

that is exactly the operator that defines canonical thermodynamical stability in the 
Newtonian Gravity [T31 [H] . We reach the rather strange and intriguing result that the 
microcanonical ensemble in General Relativity transforms to the canonical ensemble 
in the Newtonian limit! We speculate this is a hint for the presence of an intrinsic 
heat bath in General Relativity. 



6. Gravothermal catastrophe in General Relativity with A 

We are interested in this section on the effect of the cosmological constant to the 
thermodynamical stability of self-gravitating gas. Our interest is mainly theoretical 
and we basically want to see qualitatively and quantitatively how various -positive 
or negative- values of the cosmological constant affect the thermodynamical stability. 
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However, for the case of negative cosmological constant the main reason for studying 
this problem is AdS/CFT and for the case of a positive cosmological constant, 
apart from theoretical there is also physical justification for studying this problem. 
Although, at the present epoch of the evolution of the Universe the cosmological 
constant is considered to be very small, in cosmological models with a time-varying 
cosmological constant (decaying vacuum) [22j [23j [24] it is assumed to have been much 
bigger in the past so that it could have affected the formation of stars during the 
evolution of the Universe. In addition, in some types of stars, e.g. boson stars, there 
is an effective cosmological constant, i.e. a term in Einstein's equations similar to the 
cosmological constant [25j [26l [27] with a value, big enough to significantly affect the 
star's configuration. Therefore, the study of stability of TOV for various values of the 
cosmological constant is of immediate interest for Astrophysics. 

This work can be considered as the completion of previous works in the Newtonian 
case [IU [TI>J [TB] . Let us recall that gravothermal catastrophe [HI [551 12H] is called the 
set of thermodynamic instabilities in the microcanonical ensemble of a, bounded in a 
spherical box, self-gravitating gas in Newtonian gravity. Recall that the Newtonian 
limit of TOV is Emden equation [30l 03] ■ It can also be derived by entropy 
extremization, so that it describes thermodynamic equilibria as well, just like TOV 
equation (see section [3} . As we have seen for a linear equation of state 

p = qpc 2 , q = const. (59) 

the dust particles limit corresponds to q — > 0. In the present study we assume this 
linear equation of state for < q < 1. We will also assume a constant chemical 
composition and a constant entropy per nucleon (d(s/n) = 0) throughout the entire 
fluid sphere. These assumptions correspond to the following physical cases [21]: (a) 
a white dwarf or a neutron star of low mass, in which the temperature is essentially 
at absolute zero (Nernst's theorem then ensures s/n — 0); (b) stars in convective 
equilibrium, like supermassive stars and self-gravitating radiation ('photon stars'). 
Another interesting case is the 'stiff' case of matter for which q = 1 where the speed 
of sound equals the speed of light (e.g. black hole gas in Banks-Fischler cosmology 

The problem of dynamical stability under the same assumptions in the 
asymptotically flat case has been studied by Chavanis [32j [33]. He uses Weinberg's 
theorems (p. 305 in [5Tj) which state that, under the above assumptions, in a 
series of equilibria, the mass as long as the baryon number have a maximum at a 
turning point of stability. He proves that non-linear dynamical stability as defined in 
Weinberg coincides with linear dynamical stability for matter with a linear equation 
of state and with constant entropy per nucleon. SWZ [9] had proven for the 
specific case of radiation that linear dynamical stability coincides with microcanonical 
thermodynamical stability. We have proven in the previous section that, in the most 
general case under no specific assumptions, linear dynamical stability coincides with 
microcanonical thermodynamical stability. So that, all cases studied by Chavanis in 
|33] are in fact gravothermal instabilities. In addition this equivalence we have proven, 
hints to a relation between the stability theorems of Weinberg with the one of Poincare 
[Ml [33]. In fact it seems that Weinberg's theorem is just a specific case of Poincare's 
theorem of linear series of equilibria. 

Here, we are mainly interested at the type of instability that is related with the 
complete absent of equilibria. That is, with the region above the mass or baryon 
number maximum in a series of equilibria. This region is characterized by a minimum 
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radius down to which equilibrium configurations do exist. The other type of instability 
is a 'weaker' one and is related with equilibrium configurations, i.e. true equilibria in 
a series of equilibria, which nevertheless are unstable, i.e. the entropy is a saddle point 
for these equilibria. However, this instability is also characterized by the same point, 
i.e. the maximum of M and N (which coincide as Weinberg has proven) at which the 
stable equilibria become unstable. This point is called 'a turning point of stability'. 
So that, practically our scope is to determine the maximum of iV for various values 
of A in a series of equilibria, i.e. equilibria for various density contrasts \og(po/pn) 
where po is the density at the center and pr at the edge of the sphere. 

Let us prove a relation that we will need in the followings. The first law of 
thermodynamics can be expressed in the form: 

Td(-) =d{^-)+pd{-) (60) 
n n n 

Under the above assumptions the term Td{s/n) can be neglected, so that we get: 

dpc 2 — — — - — dn (61) 

n 

Integrating the above equation, taking into account the equation of state (|59p. we get 
the polytropic equation: 

p = Kn q+1 (62) 

Combining the polytropic equation (1621) with the equation of state (f59|) we get: 

Now, let our fluid sphere be bounded by perfectly reflecting and non-insulating walls. 
Let us use the dimensionless variables of Chandrasekhar IMI 



p = p e v , x = rJAnGpfr-^- (64) 
V qc z 

where y, x correspond to ip 7 £ of Chandrasekhar (we prefered to keep the notation 
relevant to the Newtonian limit [221 H] ) and po is the mass density at the center of 
the fluid sphere. Let us introduce also: 

2PA 1 ( q + V * 

where R is the radius of the sphere. Using these dimensionless variables, TOV equation 
becomes: 



A = -^, n(x) = (^Gpo^V) m(r), M = m(R) (65) 



-r- = [ ~o + 1 xe t: x \ 1 — 2 -x (66) 

^ = x 2 e-^ (67) 



with initial conditions 

1/(0) = /i(0) = 
We introduce the dimensionless energy: 



2GAf 2^ q 
Rc z q + 1 
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Let z denote the value of x at R. Integrating TOV is straightforward to calculate 
H at z and after substituting it in (|69p we get: 



2q 



zy'(z) - qz 2 e-M - f z 2 (jfez^z) + l) 



0( ^^fT 2^(z) + l (70) 



where denotes |^ 



We use the dimensionless baryon number introduced by 



Chavanis [35] which we call B: 

£? = — , N* = 4ttR ( (71) 

Using equation (1251) for iV we find in dimensionless variables: 

B(z) — 3^jrr [ Z x-e-tk ( l _ 2-JL.H - V 1 (72) 



z ^pr 7 V g+l£ 3g+l / 

For a series of equilibria, where each one is specified by the dimensionless quantity z 
that determines the radius R of the sphere and/or its central density po, the quantity 
B(z) expresses the baryon number for a fixed radius, Q{z) expresses the total mass 
for a fixed radius, while 1/Q(z) expresses the radius for a fixed mass. Instead of z, 
it is equivalent to use the so called density contrast log(po/Pfl), as a free variable to 
define the equilibria. 

In the asymptotically flat case it is trivial to generate a series of equilibria, since 
the system of equations (|fJrJ|) . (p7|) can be regarded with respect to z instead of x with 
no problem. However, in the presence of A it is false to consider equation (|66p with 
respect to z. There appears the variable A in the equations which cannot be regarded as 
constant, because it contains po- So that, in fact there are two independent variables 
(A, z). To overcome this difficulty we developed a strategy similar to the one that 
has been used in order to solve the corresponding Newtonian problem [21 [15] . The 
cosmological constant introduces a mass scale to the system M\ — pa^kR 3 - We define 
the dimensionless mass 

m ~2M A 8tt PA i?3 " Az3 {l6) 

We solve numerically the system ([66 ]) - ([67]) for a z range of values choosing this value 
for A that keeps constant the dimensionless mass rh. This operation is performed by 
a computer program we developed. In this way a consistent series of equilibria can be 
generated that corresponds to a fixed M or R or p\ according to equation ([73]) . Then, 
by generating series for various rh we can see how various critical quantities vary with 
A. 

In figure[T]we can see how, the baryon number changes with respect to the density 
contrast for some fixed values of cosmological constant in a series of equilibria. The 
fixed positive value (corresponding to rh = 1) is denoted as 'dS' and the fixed negative 
value (corresponding to rh = — 1) is denoted as 'AdS', while the case A = is denoted 
as 'Flat'. The series of equilibria are plotted for the cases of q = 1/3 (radiation or 
neutron stars) and q = 1 (stiff matter). We see that, as the cosmological constant 
increases, the system is stabilized with respect to the 'strong' instability, i.e. the one 
corresponding to no equilibria, since the turning point goes to higher N values. The 
'weak' instability, i.e. the one corresponding to unstable equilibria, occurs at lower 
density contrast values as the cosmological constant increases. 
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(a) q = 1/3 




Figure 1. The baryon number N with respect to the density contrast log(po / Pr) 
for radiation (q = 1/3) and for stiff matter (q = 1) for asymptotically flat, de Sitter 
and anti-de Sitter cases. The points C; are turning points of stability, i.e. at their 
left side the equilibria arc stable, while at their right side, unstable. At each case 
above d there are no equilibria at all. 



We would like to plot the critical radius, i.e. the minimum radius of the fluid 
sphere down to which equilibria do exist, with respect to the cosmological constant. 
However, there is the following complication in this procedure. In an asymptotically 
de Sitter space, every fluid sphere of mass M and radius R has two characteristic 
length scales: the black hole radius Rbh, that would correspond to a black hole 
horizon if R < Rbh and the cosmological radius Rc, that would correspond to a 
cosmological horizon if R > Rc (see figure HJ) . By black hole radius we do not mean 
the Schwartzschild radius GM/c 2 but the larger 'Schwartzscild-de Sitter' radius that 
corresponds to the smaller real root of the polynomial with respect to Rh- 

2GM 8nG „ 2 „ 
1 - R^ - l*"** = ° (74) 
The cosmological radius, that we denote Rc, is the larger real root of the polynomial 
(|74l) . We expect and indeed verify numerically that the critical radius R cr is always 
larger than Rbh- Thus, the system becomes unstable always before it reaches its black 
hole horizon radius and become a black hole. On the contrary, as the cosmological 
constant is increasing, i.e. Rc is decreasing, one would expect that R cr might become 
larger than Rc ■ This would mean that the sphere with R = Rc = R cr is the maximum 
sphere for which an instability can be observed, because for any bigger sphere the 
critical radius would lie beyond the cosmological horizon. 

We can determine the dimensionless mass m at which R — Rc with respect to 
Q = 2GM/Rc 2 as follows. It is easy to verify that the polynomial ([71|) . for a sphere 
of radius R, and for Rh = Rc can be written in the form: 

We see that R = Rc if: 
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Figure 2. The critical radius R cr with respect to the cosmological constant pa 
in Schwartzschild units R s = 2GM/c 2 , p s = 3M/4wR^ for a fixed mass M. The 
spheres of radius R, < i?, cr are strongly thermodynamically unstable, i.e. there 
exist no equilibria at all. For R > R cr there exist equilibria which are stable or 
unstable depending on the density contrast po/Pi?- In the same plot, are drawn 
the black hole radius Rbh and the cosmological radius Rc, as well. At points A, 
the critical radius hits the cosmological horizon. 



We verify numerically that this condition holds when the critical radius 'hits' the 
cosmological horizon. Physically R cr > Rc is meaningless and this is manifested 
mathematically (and numerically) as inability to calculate various physical quantities 
(for example N becomes complex and the computer program breaks down). So that 
a fluid sphere with R cr > Rc can be considered to be stable, since matter beyond the 
cosmological horizon cannot interact with matter inside it. In figure is plotted the 
critical radius in Schwartzschild units, together with the black hole and cosmological 
radii, with respect to the cosmological constant for a fixed mass. The critical radius 
hits the cosmological horizon at point A. Spheres with radii smaller than the critical 
radius are microcanonical thermodynamically unstable, in the sense that there are no 
thermodynamic equilibria in the microcanonical ensemble. 

Recall that the Newtonian limit of TOV equation is taken for q = kT/m p c 2 <C 1. 
Let us calculate B = N/N* for dust matter. Substituting c 2 = l/qm p (3 into equation 
(jTTj) . then taking q — > 0, equation (jTTj) becomes: 

GM/3 

This is exactly the Newtonian dimensionless temperature as defined in JT5J [T3] for 
m p = 1. Since B controls microcanonical thermodynamical stability in relativity and 
canonical thermodynamical stability in Newtonian limit |15l 113] . we see once more 
that the Newtonian limit of the relativistic microcanonical ensemble is the canonical 
ensembl^l, verifying the result of section [5] This conclusion can also be intuitively 
justified, since the canonical ensemble in Newtonian Gravity behaves qualitatively 
similar to a General Relativistic system. That is because in the Newtonian canonical 
ensemble, for a fixed mass, the thermodynamical instability is triggered as the system 

§ In addition, we verified numerically that for q = 0, i.e. non-relativistic matter, the critical value of 
z for any value of the cosmological constant equals the critical z of the Newtonian canonical ensemble 
(and not of the Newtonian microcanonical ensemble). 
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is contracted beyond a minimum radius (i.e. for small radii) [15) . just like in General 
Relativity, while in the Newtonian microcanonical ensemble, the instability sets in 
when the system is expanded beyond a maximum radius (i.e. for big radii). 

7. Conclusions 

Probably, the most interesting result of this work is the fact that for a static spherically 
symmetric perfect fluid, the Einstein's equations are completely equivalent with 
the thermodynamics of self-gravitating gas in the microcanonical ensemble. This 
equivalence does not apply only on the solutions of Einstein's equations, i.e. the 
equilibria, but also on the stability of these equilibria! This result gives further support 
on the deep connection between gravity and thermodynamics. 

A rather intriguing result we found, is that the microcanonical ensemble in 
General Relativity becomes the canonical ensemble in the Newtonian limit (dust 
matter) as we have analytically calculated in section [5] It seems, as though in General 
Relativity, there is an implicitly constructed 'heat bath' so that a microcanonical 
ensemble would look like a canonical ensemble when non-relativistic matter is 
considered! 

Finally, we have studied how the stability properties are influenced by the presence 
of a cosmological constant. The main results, we could say that are summarized in the 
characteristic figures [T] and [2] For the cases of radiation and neutron cores, the critical 
radius less than which there is no thermodynamical equilibria, is decreasing with 
increasing cosmological constant while for stiff matter it is increasing. This critical 
radius is always bigger than the black hole radius so that the system becomes unstable 
long before it reaches its black hole radius for any value of the cosmological constant. 
However, at some point the critical radius hits the cosmological horizon. Any sphere 
with bigger radius can be considered stable since matter outside the cosmological 
horizon cannot interact with matter inside. Moreover, an increase in the cosmological 
constant tends to stabilize the system thermodynamically, rendering de Sitter space, 
at least at the purely classical level, more stable than anti-de Sitter. 
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Appendix A. 

We review here Yabushita's [20] calculation of the equation that determines linear 
stability, including a cosmological constant. Using the metric (TT]), the energy 
momentum tensor ([2]) has the following components (p. 251 in Tolman [17] ) 




(A.1) 



(A.2) 
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where prime denotes differentiation with respect to r and dot with respect to i. The 
four velocity has components m m = u 1 , 0, 0). At the equilibrium u 1 = 0, the energy 
momentum tensor has only the followings non-zero components 

T ° = p , Tl = T% = Tl = -p (A.6) 

and A, v are independent of t. Therefore equations (|A.1[) - (|A.3[) become 

^=~(Pe+/W?We (A.7) 

8vrG -\ fK 1 \ 1 



,2 P^^\-f--2)+-2-^ (A.8) 

8ttG x /z/ 1 A 1 

— Pe = e" A * - + - - - + A A.9 
c y r r J r 

where the suffix e denotes quantities at the equilibrium. Equation (|A.7[) is the 
relativistic analogue of Newtonian expression for the hydrostatic equilibrium 

f = -4 ^ 

dr dr 

It is derived by equating (|A.2[) and (|A.3[) . Applying the transformation 

e -e =1 _!^!l_A r2 (A . n) 

r 3 

to equation (|A.9[) . solving with respect to v' e and substituting v' e into equation (|A.7[) ■ 
we get TOV equation: 

/</ --( 1 --\-,>,)['^— + lr i a'-r- -<-r)(] . r^_L . If,* ) ' (A.12) 



c 2 3 / \ re 2 3 

that is the relativistic equation of hydrostatic equilibria. The remaining equation 
(|A.8[) gives 

m = Airp e r 2 (A.13) 

which means that m is the mass contained inside the radius r. Equations (IA.12[) and 
(|A.13j) together with an equation of state p e = p e (p e ) fully describe the equilibria, 
provided some initial conditions are given. 

Now, let us consider small perturbations about an equilibrium configuration 

A = A e + SX , v = u e + 5v , p = p e + Sp , p = p e + Sp , u 1 = Su 1 , u° = u° e + Su° (A.14) 

with all variations depending on both r and t. To first order, equations (|A.1[) - (|A.3|) 
yield 

\e-"'8\ = ^ Up' + ^(Sp + Spc 2 ) + ~(p e + Pe c 2 )sA (A.15) 
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The conservation of mass 2q = gives 



n 1 + 1 -kp+ P j)+n( l -\> + l - V ' 



(A.16) 



which after replacing (|A.14|) and differentiation with respect to t and use of equations 
(|A~7l) -(TA~9t gives 



dp 



d 2 \ < 

57rG V 9r r . 



-SX = 



,l/ e -Ae 



OP + ^(0P + OpC) + —{Pe 



PeC 



Substituting equation (jA.15[) into (|A.17[) we get 

6 -(— + - 

\dr r 

Let 

f(r,t) ee 5m(r,t) 

with 

s 1 g/ 

47rr 2 <9r 
Assuming a perturbation 

s at 

dm ~ e 

with 5m(r = 0) = and integrating equation (|A.15|) we get finally 



& + —{5p + 5pc 2 ) + —{ Pe + Pe c 2 ) = — T 



e 2 f 



(A.17) 

(A.18) 
(A.19) 
(A.20) 
(A.21) 

(A.22) 



The cosmological constant does not enter explicitly in the equation but only implicitly 
in the equilibrium quantities. 
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